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Summary. The effects of a deep-burrowing earthworm species, Martiodrilus carimaguen- 
sis, on soil fertility and plant growth were assessed in two grasslands in the Eastern Plains of 
Colombia. This species aestivates deep in the soil during the dry season. Fibre glass netting 
was used to prevent colonisation of the top soil at the onset of the humid season. Effects of 
earthworm absence on soil properties and plant biomass. were analysed with principal com- 
ponent analysis. 

Experimental manipulation succeeded in preventing M. carimaguensis activity in the top 
soil. The larger part of the variance in our data (29.35 %) was explained by the effect of pa- 
sture type and age on soil structure and organic matter, A lesser part (10.81%) was due to the 
effect of legumes on soil physical properties. The effect of earthworm absence on soil pro- 
perties and plant growth explained 16.50 % of the total variance. M. carimaguensis biomass 
was associated with low soil compaction, high C contents, low Al saturation, high herbaceous 
biomass and low weed biomass. These results support the general knowledge of how earth- 
worms can affect soil fertility and plant growth, f 
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Introduction 


Most studies of the impact of earthworms on soil have been conducted in short-term pot ex- 
periments (e.g, Pashanasi et al. 1992; Spain et al. 1992; Derouard et al. 1997). Despite the sig- 
nificant advances this approach has made in the understanding of the role of earthworms in 
soil processes, there is still a lack of understanding of how these organisms may influence 
ecosystem processes in in-situ field conditions (Bohlen et al. 1995). 

Our objective was to assess an innovative method (Baker 1998) to selectively exclude ane- 
cic earthworms (sensu Bouché 1977) from small experimental units. The efficiency of the 
method is discussed in relation to its ability to eliminate the large anecic earthworm Marti- 
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odrilus carimaguensis Jiménez and Moreno from the top-soil of two grass-legume based 
pastures in the Eastern Plains of Colombia. At Carimagua, populations of this species are en- 
hanced in intensive pastures and dramatically depleted in annual crops (Jiménez et al. 1998a). 
The effects of the disappearance of this species on soil properties and plant growth was as- 
sessed and compared with the effects of other environmental factors. 


Materials and Methods 


Study site 


The study was carried out at the CIAT-CORPOICA research station of Carimagua, Eastern Plains, Meta, 
Colombia (4° 37'N, 71° 19'W). The climate is humid tropical with an annual mean temperature and 
rainfall of 26°C and 2300 mm respectively (CIAT data). Native vegetation is a patchwork of open 
savannas and gallery forests. Soils are well aggregated Oxisols with high acidity and very low chemical 
fertility. 


Experimental plots and design 


This experiment was carried out in two pastures cach grazed by cattle with 2.0 animal unit. ha-t: (A) A 
3 year-old pasture of Brachiaria humidicola, Arachis pintoi, Stylosauthes capitata and Centrocema acu- 
tifolium; (B) A one year-old pasture of Panicum maximum and Arachis pintoi. 

The background ecology of the earthworm species of Carimagua are available in Jiménez et al. 
(1998a, b). All species aestivate during the dry season. The anecic M. carimaguensis burrows deep into 
the soil profile (60-100 cm) and falls into diapause until the next wet season (Jiménez et al. 1998b). Our 
idea was to prevent individuals from returning to the top soil at the onset of the rains. As other species 
aestivate in superficial soil layers, we expected this experiment to have specific effects on M. carima- 
guensis. 

Twelve small experimental units were established in each pasture at the end of the dry season (early 
March 1996) and submitted to grazing pressure. The protocol is synthesised in Fig. 1: (1) a 40 x 40 x 
30 cm soil monolith was dug out with a spade; (2) the monolith sides (except the top side) were wrap- 
ped up in a fibre glass netting of 0.5 mm mesh; (3) the monolith was placed in a hole, and surface 
earthworm movements were avoided by surrounding the units with a buried metallic plate; (4) an area 
was delimited where individuals of M. carimaguensis (located by the presence of fresh casts, Jiménez 
et al. 1998b) were systematically controlled by injecting 0.01 g of active carbofuran in their burrow. 
Twelve plots were also established without fibre glass netting as controls. Plot location was randomly 
chosen in the field, avoiding the field edges. 

Soil and vegetation were sampled before, 6 months and 18 months after the experiment was 
established (respectively October 1995, August 1996 and June 1997). At each date, four units were ran- 
domly chosen in each pasture. Samples aimed at characterising vegetation (shoot and root biomass), soil 
properties (bulk density, aggregate mean weigh density, penetration resistance, sorptivity, infiltration, 
chemical properties) and earthworm (density and biomass), ant and termite populations (number of 
chambers). Methods were standardised methods recommended by the TSBF (Tropical Soil Biology and 
Fertility) program (Anderson & Ingram 1989). All samples’ were repeated 4 fold. 


Statistics 


The effect of the experimental manipulation on both earthworm biomass and the abundance of ant and 
termite chambers was tested using the Fisher PLSD test. Asymmetry of data frequency distribution was 
reduced by Box-Cox transformation (Sokal & Rohlf 1995) with the software Vernorm (Legendre & 
Vaudor 1991). The Bonferonni procedure for nested tests was used and the adjusted 0.05, 0.01 and 0.001 
significance levels were respectively: 0.008, 0.002 and 0.0002. 

A Principal Component Analysis (PCA) was performed with ADE software (Thioulouse et al. 1997) 
to extract the principal source of variance in the environmental data sets. We used a matrix of 40 lines 
(samples) x 23 columns (soil and vegetation variables). The biomass of each earthworm species and the 
abundance of ant and termite chambers were projected on the PCA axes as additional columns. Simple 
linear regressions were used to test the correlations between faunal parameters and PCA axes. : 
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Fig. 1. Description of the sampling procedure. Numbers correspond to the chronological order of the 
different manipulation steps, and refer to the full descriptions in the text 


Results 


Only one individual of M. carimaguensis was found in the experimental units at the end of the 
study. Due to the large size of M. carimaguensis, slight decreases in the density of this spe- 
cies resulted in sharp decreases in overall earthworm biomass (from ~ 89% to- 100%,p < 
0.01) (Fig. 2a, 2b). Differences observed for other earthworm species (Fig. 2a, 2b) and social 
arthropods (not shown) were not significant. 

Axis 1 of the PCA accounted for 29.33 % of the total inertia (Table 1). The correlation ratio 
associated with the variables showed negative relationships between (i) high graminea, herb- 
aceous and root biomass, high bulk density in the 5-10 cm layer, high penetration resistance, 
high pH and available nutrient content, and (ii) high aluminium saturation (Fig. 3a, Table 1). 
Ordination of the sample scores showed a separation between (i) samples taken in the first 
two dates and those taken later (Fig. 3b), and (ii) samples taken in the pasture A and in the 
pasture B (Fig: 3c). This axis was interpreted as the effect of pasture characteristics and age 
on soil physical and chemical properties. 

Axis 2 accounted for 16.52 % of the total inertia (Table 1). It opposed (i) high weed bio- 
mass, high bulk density, high penetration resistance, high total N content and high aluminium 
Saturation to (ii) high herbaceous (legume and graminea) and root biomass, high total C con- 
tent and high mean diameter of soil aggregates (Fig. 3a, Table 1). It clearly segregated 
samples taken inside the exclusion units with those taken in the control units (Fig. 3d), and 
was defined as the effect of the exclusion of M. carimaguensis on soil fertility and plant 
growth. 

Axis 3 accounted for 10.81 % of the total inertia. It was driven by the abundance of legu- 
mes in the plant cover and the soil hydraulic properties (sorptivity and infiltration) (Table 1). 
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Fig. 2. Density (a) and biomass (b) 
of each earthworm species in the 
experimental units with (+) and 
without (—) netting, at each sam- 
pling date and in the two pastures 
(N = 8 for Om and 4 for the other 
sampling dates; ind. = individuals; 
fm = fresh matter; Om = before; 
6 m= 6 months; 18 m = 18 months) 


Table 1. Relative inertia for each of the three first axes (in brackets) and correlation ratios associated 
with the variables included in the principal correspondence analysis (PCA) 





Variables Code Axis | Axis 2 Axis 3 
(29.33 %) (16.52%) (10.81 %) 
Vegetation characteristics 
Graminea biomass GB -0.751 —0.447 -0.039 
Legume biomass LB -0.173 —-0.402 -0.709 
Weed biomass WB -0.086 0.492 -0.044 
% of legumes HL 0.080 -0.101 -0.796 
Herbaceous biomass HB —0.748 -0,464 -0.191 
Root biomass RB -0.707 —0.500 -0.245 
Soil chemical properties (0-5 cm) 
pH pH -0.786 0.159 0.114 
Total C content C% —0.168 -0.359 0.018 
Total N content N% ~0.188 0.532 -0.089 
Total P content TP —-0.738 0.190 0.030 
Al saturation Al 0.625 0.493 -0.360 
Ca content Ca -0.817 -0.332 0.217 
Mg content Mg -0.803 —0.347 0.262 
K content K —0.336 —0.182 0.251 
Soil physical properties 
Bulk density (0-5 cm) BDO-5 -0.281 0.518 -0.131 
Bulk density (5-10 cm) BD5-10 -0.612 0.515 -0.091 
Mean weigh diameter of aggr. (0—10 cm) MWD -0.054 —0.377 0.167 
Penetration resistance (0-15 cm) PR0-15 —0.623 0.538 0.205 
Penetration resistance (15-30 cm) PR15-30 -0.651 0.499 0.181 
Penetration resistance (30-60 cm) PR30-60 ~0.464 0.400 0.198 
Penetration resistance (0-60 cm) PRO-60 -0.703 0.579 —0.021 
Sorptivity 5 0.164 —0.171 -0.645 
Infiltration constant A -0.174 0.030 -0.579 
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Fig. 3. Results of the principal component analysis (PCA): (a) correlation circle and diagram of axis 
eigenvalues (variables in light letters; additional variables in bold letters; codes as in Table 1; Mc = 
Martiodrilus carimaguensis biomass; Ay = Andiodrilus yoparensis biomass; Gs = Glossodrilus stkuani 
biomass; Oc = Ocnerodrilidae sp. biomass; Ant = ant abundance; Term = termite abundance); (b) (c) 
and (d) ordination of the samples in the plane defined by the axis 1 and 2 of the PCA (insert in the 
fig. 3b represents the scales of the fig. 3b, 3c and 3d; 0m = before; 6m = 6 months; 18m = 18 months; 
A= B. humidicola pasture, B = P. maximum pasture; + = with netting; — = without netting) 


This axis was interpreted as the effect of legumes on soil physical properties. Samples were 
not ordinated along this axis according to a logical scheme as was the case for the first two 
axes (not shown). 

When projected on the plan defined by the two first axes of the PCA, the biomass of M. ca- 
rimaguensis was negatively correlated with axis 2 (Fig. 3a). This negative relationship was 
confirmed by the regression coefficients (r = 0.63; p < 0.001). The biomass of other earth 
worm species and the abundance of social insect chambers were located near the axes inter- 
sections, showing that these variables were not correlated with any of them (Fig. 3a). No sig- 
nificant correlations were found between faunal biomass and the third PCA axis. 
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Discussion 


Field experiments are pertinent for evaluating earthworm effects on soil parameters and plant 
growth. This may be very useful in assessing the role of earthworms on primary production 
in both natural and managed ecosystems (Bohlen et al. 1995). In our study, preventing verti- 
cal movements of earthworm at an appropriate time and depth was an efficient method to eli- 
minate a targeted anecic species. 

The main factors influencing soil properties were pasture type and age. This may be due to 
plant specific effects, soil protection by the plant cover, litter quality and production, and the 
cumulative effects of animal trampling that may be different according to the history and the 
nature of pastures (Sanchez et al. 1991; Gijsman 1996; Gijsman & Thomas 1996). The effect 
of legumes on soil physical properties was of third importance. This may be explained by the 
effect of root channels that increase soil macroporosity and permeability (Mytton et al. 1993). 

Disappearance of M. carimaguensis resulted in significant changes in soil properties and 
plant biomass, although this impact was of secondary importance when compared with 
pasture effect. The decrease in earthworm biomass was associated with soil degradation (in- 
creased soil compaction and aluminium saturation, decreased carbon content) and herbaceous 
biomass. Decrease of plant biomass, in turn, favoured opportunistic weed species. These re- 
sults are consistent with the general knowledge of how earthworms can influence soil pro- 
perties and plant growth (see reviews by Lee 1985; Lal 1988; Edwards & Bohlen 1996; La- 
velle 1997). It demonstrates that the loss of one species in an earthworm community, when it 
is associated with an important decrease in biomass, may result in significant losses in eco- 
system function. The rapidity and intensity of the effects observed in our study, however, may 
have been accentuated by intensive animal trampling. in natural undisturbed conditions, 
earthworm impacts on soil processes are known to be persistent on a larger time scale (i.e. se- 
veral months to years) (Blanchart 1992). 

Earthworm populations are highly sensitive to land use practices (see, e.g. Lavelle & Pas- 
hansi 1989; Reddy et al. 1995), At Carimagua, native species biomass is greatly increased in 
man-made pasture or, in contrast, dramatically decreased in annual crops (Decaéns et al. 
1994; Jiménez ct al, 1998a). Such impacts may have determinant consequences on soil func- 
tion and primary production (Lavelle 1996; Lavelle 1997). Attention should now be paid to 
managing earthworm populations in tropical agroecosystems in order to profit from their im- 
pacts on soil fertility and enhance the sustainability of agricultural production. 
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